\(\int \frac {1}{\sqrt {d+e x} (a+c x^2)^3} \, dx\) [647]

Optimal result

Integrand size = 19, antiderivative size = 920 \[ \int \frac {1}{\sqrt {d+e x} \left (a+c x^2\right )^3} \, dx=\frac {(a e+c d x) \sqrt {d+e x}}{4 a \left (c d^2+a e^2\right ) \left (a+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (a e \left (c d^2+7 a e^2\right )+6 c d \left (c d^2+2 a e^2\right ) x\right )}{16 a^2 \left (c d^2+a e^2\right )^2 \left (a+c x^2\right )}+\frac {3 e \left (2 c^2 d^4+5 a c d^2 e^2+7 a^2 e^4+2 \sqrt {c} d \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )\right ) \text {arctanh}\left (\frac {\sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}-\sqrt {2} \sqrt [4]{c} \sqrt {d+e x}}{\sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}\right )}{32 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}-\frac {3 e \left (2 c^2 d^4+5 a c d^2 e^2+7 a^2 e^4+2 \sqrt {c} d \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )\right ) \text {arctanh}\left (\frac {\sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}+\sqrt {2} \sqrt [4]{c} \sqrt {d+e x}}{\sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}\right )}{32 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}-\frac {3 e \left (2 c^2 d^4+5 a c d^2 e^2+7 a^2 e^4-2 \sqrt {c} d \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )\right ) \log \left (\sqrt {c d^2+a e^2}-\sqrt {2} \sqrt [4]{c} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \sqrt {d+e x}+\sqrt {c} (d+e x)\right )}{64 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}+\frac {3 e \left (2 c^2 d^4+5 a c d^2 e^2+7 a^2 e^4-2 \sqrt {c} d \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )\right ) \log \left (\sqrt {c d^2+a e^2}+\sqrt {2} \sqrt [4]{c} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \sqrt {d+e x}+\sqrt {c} (d+e x)\right )}{64 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}} \]

[Out]

1/4*(c*d*x+a*e)*(e*x+d)^(1/2)/a/(a*e^2+c*d^2)/(c*x^2+a)^2+1/16*(a*e*(7*a*e^2+c*d^2)+6*c*d*(2*a*e^2+c*d^2)*x)*(
e*x+d)^(1/2)/a^2/(a*e^2+c*d^2)^2/(c*x^2+a)+3/64*e*arctanh((-c^(1/4)*2^(1/2)*(e*x+d)^(1/2)+(d*c^(1/2)+(a*e^2+c*
d^2)^(1/2))^(1/2))/(d*c^(1/2)-(a*e^2+c*d^2)^(1/2))^(1/2))*(2*c^2*d^4+5*a*c*d^2*e^2+7*a^2*e^4+2*d*(2*a*e^2+c*d^
2)*c^(1/2)*(a*e^2+c*d^2)^(1/2))/a^2/c^(1/4)/(a*e^2+c*d^2)^(5/2)*2^(1/2)/(d*c^(1/2)-(a*e^2+c*d^2)^(1/2))^(1/2)-
3/64*e*arctanh((c^(1/4)*2^(1/2)*(e*x+d)^(1/2)+(d*c^(1/2)+(a*e^2+c*d^2)^(1/2))^(1/2))/(d*c^(1/2)-(a*e^2+c*d^2)^
(1/2))^(1/2))*(2*c^2*d^4+5*a*c*d^2*e^2+7*a^2*e^4+2*d*(2*a*e^2+c*d^2)*c^(1/2)*(a*e^2+c*d^2)^(1/2))/a^2/c^(1/4)/
(a*e^2+c*d^2)^(5/2)*2^(1/2)/(d*c^(1/2)-(a*e^2+c*d^2)^(1/2))^(1/2)-3/128*e*ln((e*x+d)*c^(1/2)+(a*e^2+c*d^2)^(1/
2)-c^(1/4)*2^(1/2)*(e*x+d)^(1/2)*(d*c^(1/2)+(a*e^2+c*d^2)^(1/2))^(1/2))*(2*c^2*d^4+5*a*c*d^2*e^2+7*a^2*e^4-2*d
*(2*a*e^2+c*d^2)*c^(1/2)*(a*e^2+c*d^2)^(1/2))/a^2/c^(1/4)/(a*e^2+c*d^2)^(5/2)*2^(1/2)/(d*c^(1/2)+(a*e^2+c*d^2)
^(1/2))^(1/2)+3/128*e*ln((e*x+d)*c^(1/2)+(a*e^2+c*d^2)^(1/2)+c^(1/4)*2^(1/2)*(e*x+d)^(1/2)*(d*c^(1/2)+(a*e^2+c
*d^2)^(1/2))^(1/2))*(2*c^2*d^4+5*a*c*d^2*e^2+7*a^2*e^4-2*d*(2*a*e^2+c*d^2)*c^(1/2)*(a*e^2+c*d^2)^(1/2))/a^2/c^
(1/4)/(a*e^2+c*d^2)^(5/2)*2^(1/2)/(d*c^(1/2)+(a*e^2+c*d^2)^(1/2))^(1/2)

Rubi [A] (verified)

Time = 4.18 (sec) , antiderivative size = 920, normalized size of antiderivative = 1.00, number of steps used = 12, number of rules used = 8, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.421, Rules used = {755, 837, 841, 1183, 648, 632, 212, 642} \[ \int \frac {1}{\sqrt {d+e x} \left (a+c x^2\right )^3} \, dx=\frac {\sqrt {d+e x} (a e+c d x)}{4 a \left (c d^2+a e^2\right ) \left (c x^2+a\right )^2}+\frac {3 e \left (2 c^2 d^4+5 a c e^2 d^2+2 \sqrt {c} \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right ) d+7 a^2 e^4\right ) \text {arctanh}\left (\frac {\sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}-\sqrt {2} \sqrt [4]{c} \sqrt {d+e x}}{\sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}\right )}{32 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}-\frac {3 e \left (2 c^2 d^4+5 a c e^2 d^2+2 \sqrt {c} \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right ) d+7 a^2 e^4\right ) \text {arctanh}\left (\frac {\sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}+\sqrt {2} \sqrt [4]{c} \sqrt {d+e x}}{\sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}\right )}{32 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}-\frac {3 e \left (2 c^2 d^4+5 a c e^2 d^2-2 \sqrt {c} \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right ) d+7 a^2 e^4\right ) \log \left (\sqrt {c} (d+e x)-\sqrt {2} \sqrt [4]{c} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \sqrt {d+e x}+\sqrt {c d^2+a e^2}\right )}{64 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}+\frac {3 e \left (2 c^2 d^4+5 a c e^2 d^2-2 \sqrt {c} \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right ) d+7 a^2 e^4\right ) \log \left (\sqrt {c} (d+e x)+\sqrt {2} \sqrt [4]{c} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \sqrt {d+e x}+\sqrt {c d^2+a e^2}\right )}{64 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}+\frac {\sqrt {d+e x} \left (a e \left (c d^2+7 a e^2\right )+6 c d \left (c d^2+2 a e^2\right ) x\right )}{16 a^2 \left (c d^2+a e^2\right )^2 \left (c x^2+a\right )} \]

[In]

Int[1/(Sqrt[d + e*x]*(a + c*x^2)^3),x]

[Out]

((a*e + c*d*x)*Sqrt[d + e*x])/(4*a*(c*d^2 + a*e^2)*(a + c*x^2)^2) + (Sqrt[d + e*x]*(a*e*(c*d^2 + 7*a*e^2) + 6*
c*d*(c*d^2 + 2*a*e^2)*x))/(16*a^2*(c*d^2 + a*e^2)^2*(a + c*x^2)) + (3*e*(2*c^2*d^4 + 5*a*c*d^2*e^2 + 7*a^2*e^4
+ 2*Sqrt[c]*d*Sqrt[c*d^2 + a*e^2]*(c*d^2 + 2*a*e^2))*ArcTanh[(Sqrt[Sqrt[c]*d + Sqrt[c*d^2 + a*e^2]] - Sqrt[2]
*c^(1/4)*Sqrt[d + e*x])/Sqrt[Sqrt[c]*d - Sqrt[c*d^2 + a*e^2]]])/(32*Sqrt[2]*a^2*c^(1/4)*(c*d^2 + a*e^2)^(5/2)*
Sqrt[Sqrt[c]*d - Sqrt[c*d^2 + a*e^2]]) - (3*e*(2*c^2*d^4 + 5*a*c*d^2*e^2 + 7*a^2*e^4 + 2*Sqrt[c]*d*Sqrt[c*d^2
+ a*e^2]*(c*d^2 + 2*a*e^2))*ArcTanh[(Sqrt[Sqrt[c]*d + Sqrt[c*d^2 + a*e^2]] + Sqrt[2]*c^(1/4)*Sqrt[d + e*x])/Sq
rt[Sqrt[c]*d - Sqrt[c*d^2 + a*e^2]]])/(32*Sqrt[2]*a^2*c^(1/4)*(c*d^2 + a*e^2)^(5/2)*Sqrt[Sqrt[c]*d - Sqrt[c*d^
2 + a*e^2]]) - (3*e*(2*c^2*d^4 + 5*a*c*d^2*e^2 + 7*a^2*e^4 - 2*Sqrt[c]*d*Sqrt[c*d^2 + a*e^2]*(c*d^2 + 2*a*e^2)
)*Log[Sqrt[c*d^2 + a*e^2] - Sqrt[2]*c^(1/4)*Sqrt[Sqrt[c]*d + Sqrt[c*d^2 + a*e^2]]*Sqrt[d + e*x] + Sqrt[c]*(d +
e*x)])/(64*Sqrt[2]*a^2*c^(1/4)*(c*d^2 + a*e^2)^(5/2)*Sqrt[Sqrt[c]*d + Sqrt[c*d^2 + a*e^2]]) + (3*e*(2*c^2*d^4
+ 5*a*c*d^2*e^2 + 7*a^2*e^4 - 2*Sqrt[c]*d*Sqrt[c*d^2 + a*e^2]*(c*d^2 + 2*a*e^2))*Log[Sqrt[c*d^2 + a*e^2] + Sq
rt[2]*c^(1/4)*Sqrt[Sqrt[c]*d + Sqrt[c*d^2 + a*e^2]]*Sqrt[d + e*x] + Sqrt[c]*(d + e*x)])/(64*Sqrt[2]*a^2*c^(1/4
)*(c*d^2 + a*e^2)^(5/2)*Sqrt[Sqrt[c]*d + Sqrt[c*d^2 + a*e^2]])

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 632

Int[((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> Dist[-2, Subst[Int[1/Simp[b^2 - 4*a*c - x^2, x], x]
, x, b + 2*c*x], x] /; FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]

Rule 642

Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Simp[d*(Log[RemoveContent[a + b*x +
c*x^2, x]]/b), x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[2*c*d - b*e, 0]

Rule 648

Int[((d_.) + (e_.)*(x_))/((a_) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Dist[(2*c*d - b*e)/(2*c), Int[1/(a +
b*x + c*x^2), x], x] + Dist[e/(2*c), Int[(b + 2*c*x)/(a + b*x + c*x^2), x], x] /; FreeQ[{a, b, c, d, e}, x] &
& NeQ[2*c*d - b*e, 0] && NeQ[b^2 - 4*a*c, 0] &&  !NiceSqrtQ[b^2 - 4*a*c]

Rule 755

Int[((d_) + (e_.)*(x_))^(m_)*((a_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(-(d + e*x)^(m + 1))*(a*e + c*d*x)*
((a + c*x^2)^(p + 1)/(2*a*(p + 1)*(c*d^2 + a*e^2))), x] + Dist[1/(2*a*(p + 1)*(c*d^2 + a*e^2)), Int[(d + e*x)^
m*Simp[c*d^2*(2*p + 3) + a*e^2*(m + 2*p + 3) + c*e*d*(m + 2*p + 4)*x, x]*(a + c*x^2)^(p + 1), x], x] /; FreeQ[
{a, c, d, e, m}, x] && NeQ[c*d^2 + a*e^2, 0] && LtQ[p, -1] && IntQuadraticQ[a, 0, c, d, e, m, p, x]

Rule 837

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(-(d + e*x)^(
m + 1))*(f*a*c*e - a*g*c*d + c*(c*d*f + a*e*g)*x)*((a + c*x^2)^(p + 1)/(2*a*c*(p + 1)*(c*d^2 + a*e^2))), x] +
Dist[1/(2*a*c*(p + 1)*(c*d^2 + a*e^2)), Int[(d + e*x)^m*(a + c*x^2)^(p + 1)*Simp[f*(c^2*d^2*(2*p + 3) + a*c*e^
2*(m + 2*p + 3)) - a*c*d*e*g*m + c*e*(c*d*f + a*e*g)*(m + 2*p + 4)*x, x], x], x] /; FreeQ[{a, c, d, e, f, g},
x] && NeQ[c*d^2 + a*e^2, 0] && LtQ[p, -1] && (IntegerQ[m] || IntegerQ[p] || IntegersQ[2*m, 2*p])

Rule 841

Int[((f_.) + (g_.)*(x_))/(Sqrt[(d_.) + (e_.)*(x_)]*((a_) + (c_.)*(x_)^2)), x_Symbol] :> Dist[2, Subst[Int[(e*f
- d*g + g*x^2)/(c*d^2 + a*e^2 - 2*c*d*x^2 + c*x^4), x], x, Sqrt[d + e*x]], x] /; FreeQ[{a, c, d, e, f, g}, x]
&& NeQ[c*d^2 + a*e^2, 0]

Rule 1183

Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[a/c, 2]}, With[{r =
Rt[2*q - b/c, 2]}, Dist[1/(2*c*q*r), Int[(d*r - (d - e*q)*x)/(q - r*x + x^2), x], x] + Dist[1/(2*c*q*r), Int[(
d*r + (d - e*q)*x)/(q + r*x + x^2), x], x]]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2
- b*d*e + a*e^2, 0] && NegQ[b^2 - 4*a*c]

Rubi steps \begin{align*} \text {integral}& = \frac {(a e+c d x) \sqrt {d+e x}}{4 a \left (c d^2+a e^2\right ) \left (a+c x^2\right )^2}-\frac {\int \frac {\frac {1}{2} \left (-6 c d^2-7 a e^2\right )-\frac {5}{2} c d e x}{\sqrt {d+e x} \left (a+c x^2\right )^2} \, dx}{4 a \left (c d^2+a e^2\right )} \\ & = \frac {(a e+c d x) \sqrt {d+e x}}{4 a \left (c d^2+a e^2\right ) \left (a+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (a e \left (c d^2+7 a e^2\right )+6 c d \left (c d^2+2 a e^2\right ) x\right )}{16 a^2 \left (c d^2+a e^2\right )^2 \left (a+c x^2\right )}+\frac {\int \frac {\frac {3}{4} c \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )+\frac {3}{2} c^2 d e \left (c d^2+2 a e^2\right ) x}{\sqrt {d+e x} \left (a+c x^2\right )} \, dx}{8 a^2 c \left (c d^2+a e^2\right )^2} \\ & = \frac {(a e+c d x) \sqrt {d+e x}}{4 a \left (c d^2+a e^2\right ) \left (a+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (a e \left (c d^2+7 a e^2\right )+6 c d \left (c d^2+2 a e^2\right ) x\right )}{16 a^2 \left (c d^2+a e^2\right )^2 \left (a+c x^2\right )}+\frac {\text {Subst}\left (\int \frac {-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )+\frac {3}{2} c^2 d e \left (c d^2+2 a e^2\right ) x^2}{c d^2+a e^2-2 c d x^2+c x^4} \, dx,x,\sqrt {d+e x}\right )}{4 a^2 c \left (c d^2+a e^2\right )^2} \\ & = \frac {(a e+c d x) \sqrt {d+e x}}{4 a \left (c d^2+a e^2\right ) \left (a+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (a e \left (c d^2+7 a e^2\right )+6 c d \left (c d^2+2 a e^2\right ) x\right )}{16 a^2 \left (c d^2+a e^2\right )^2 \left (a+c x^2\right )}+\frac {\text {Subst}\left (\int \frac {\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \left (-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right )}{\sqrt [4]{c}}-\left (-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )-\frac {3}{2} c^{3/2} d e \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right ) x}{\frac {\sqrt {c d^2+a e^2}}{\sqrt {c}}-\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} x}{\sqrt [4]{c}}+x^2} \, dx,x,\sqrt {d+e x}\right )}{8 \sqrt {2} a^2 c^{5/4} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}+\frac {\text {Subst}\left (\int \frac {\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \left (-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right )}{\sqrt [4]{c}}+\left (-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )-\frac {3}{2} c^{3/2} d e \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right ) x}{\frac {\sqrt {c d^2+a e^2}}{\sqrt {c}}+\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} x}{\sqrt [4]{c}}+x^2} \, dx,x,\sqrt {d+e x}\right )}{8 \sqrt {2} a^2 c^{5/4} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}} \\ & = \frac {(a e+c d x) \sqrt {d+e x}}{4 a \left (c d^2+a e^2\right ) \left (a+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (a e \left (c d^2+7 a e^2\right )+6 c d \left (c d^2+2 a e^2\right ) x\right )}{16 a^2 \left (c d^2+a e^2\right )^2 \left (a+c x^2\right )}+\frac {\left (\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )+\frac {3}{2} c^{3/2} d e \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )-\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right ) \text {Subst}\left (\int \frac {-\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}{\sqrt [4]{c}}+2 x}{\frac {\sqrt {c d^2+a e^2}}{\sqrt {c}}-\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} x}{\sqrt [4]{c}}+x^2} \, dx,x,\sqrt {d+e x}\right )}{16 \sqrt {2} a^2 c^{5/4} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}+\frac {\left (-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )-\frac {3}{2} c^{3/2} d e \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right ) \text {Subst}\left (\int \frac {\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}{\sqrt [4]{c}}+2 x}{\frac {\sqrt {c d^2+a e^2}}{\sqrt {c}}+\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} x}{\sqrt [4]{c}}+x^2} \, dx,x,\sqrt {d+e x}\right )}{16 \sqrt {2} a^2 c^{5/4} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}+\frac {\left (\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \left (\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )+\frac {3}{2} c^{3/2} d e \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )-\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right )}{\sqrt [4]{c}}+\frac {2 \sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \left (-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right )}{\sqrt [4]{c}}\right ) \text {Subst}\left (\int \frac {1}{\frac {\sqrt {c d^2+a e^2}}{\sqrt {c}}-\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} x}{\sqrt [4]{c}}+x^2} \, dx,x,\sqrt {d+e x}\right )}{16 \sqrt {2} a^2 c^{5/4} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}+\frac {\left (\frac {2 \sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \left (-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right )}{\sqrt [4]{c}}-\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \left (-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )-\frac {3}{2} c^{3/2} d e \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right )}{\sqrt [4]{c}}\right ) \text {Subst}\left (\int \frac {1}{\frac {\sqrt {c d^2+a e^2}}{\sqrt {c}}+\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} x}{\sqrt [4]{c}}+x^2} \, dx,x,\sqrt {d+e x}\right )}{16 \sqrt {2} a^2 c^{5/4} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}} \\ & = \frac {(a e+c d x) \sqrt {d+e x}}{4 a \left (c d^2+a e^2\right ) \left (a+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (a e \left (c d^2+7 a e^2\right )+6 c d \left (c d^2+2 a e^2\right ) x\right )}{16 a^2 \left (c d^2+a e^2\right )^2 \left (a+c x^2\right )}-\frac {3 e \left (2 c^2 d^4+5 a c d^2 e^2+7 a^2 e^4-\sqrt {c} d \sqrt {c d^2+a e^2} \left (2 c d^2+4 a e^2\right )\right ) \log \left (\sqrt {c d^2+a e^2}-\sqrt {2} \sqrt [4]{c} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \sqrt {d+e x}+\sqrt {c} (d+e x)\right )}{64 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}+\frac {3 e \left (2 c^2 d^4+5 a c d^2 e^2+7 a^2 e^4-\sqrt {c} d \sqrt {c d^2+a e^2} \left (2 c d^2+4 a e^2\right )\right ) \log \left (\sqrt {c d^2+a e^2}+\sqrt {2} \sqrt [4]{c} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \sqrt {d+e x}+\sqrt {c} (d+e x)\right )}{64 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}-\frac {\left (\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \left (\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )+\frac {3}{2} c^{3/2} d e \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )-\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right )}{\sqrt [4]{c}}+\frac {2 \sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \left (-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right )}{\sqrt [4]{c}}\right ) \text {Subst}\left (\int \frac {1}{2 \left (d-\frac {\sqrt {c d^2+a e^2}}{\sqrt {c}}\right )-x^2} \, dx,x,-\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}{\sqrt [4]{c}}+2 \sqrt {d+e x}\right )}{8 \sqrt {2} a^2 c^{5/4} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}-\frac {\left (\frac {2 \sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \left (-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right )}{\sqrt [4]{c}}-\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \left (-\frac {3}{2} c^2 d^2 e \left (c d^2+2 a e^2\right )-\frac {3}{2} c^{3/2} d e \sqrt {c d^2+a e^2} \left (c d^2+2 a e^2\right )+\frac {3}{4} c e \left (4 c^2 d^4+9 a c d^2 e^2+7 a^2 e^4\right )\right )}{\sqrt [4]{c}}\right ) \text {Subst}\left (\int \frac {1}{2 \left (d-\frac {\sqrt {c d^2+a e^2}}{\sqrt {c}}\right )-x^2} \, dx,x,\frac {\sqrt {2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}{\sqrt [4]{c}}+2 \sqrt {d+e x}\right )}{8 \sqrt {2} a^2 c^{5/4} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}} \\ & = \frac {(a e+c d x) \sqrt {d+e x}}{4 a \left (c d^2+a e^2\right ) \left (a+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (a e \left (c d^2+7 a e^2\right )+6 c d \left (c d^2+2 a e^2\right ) x\right )}{16 a^2 \left (c d^2+a e^2\right )^2 \left (a+c x^2\right )}+\frac {3 e \left (2 c^2 d^4+5 a c d^2 e^2+7 a^2 e^4+\sqrt {c} d \sqrt {c d^2+a e^2} \left (2 c d^2+4 a e^2\right )\right ) \tanh ^{-1}\left (\frac {\sqrt [4]{c} \left (\frac {\sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}{\sqrt [4]{c}}-\sqrt {2} \sqrt {d+e x}\right )}{\sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}\right )}{32 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}-\frac {3 e \left (2 c^2 d^4+5 a c d^2 e^2+7 a^2 e^4+\sqrt {c} d \sqrt {c d^2+a e^2} \left (2 c d^2+4 a e^2\right )\right ) \tanh ^{-1}\left (\frac {\sqrt [4]{c} \left (\frac {\sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}{\sqrt [4]{c}}+\sqrt {2} \sqrt {d+e x}\right )}{\sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}\right )}{32 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d-\sqrt {c d^2+a e^2}}}-\frac {3 e \left (2 c^2 d^4+5 a c d^2 e^2+7 a^2 e^4-\sqrt {c} d \sqrt {c d^2+a e^2} \left (2 c d^2+4 a e^2\right )\right ) \log \left (\sqrt {c d^2+a e^2}-\sqrt {2} \sqrt [4]{c} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \sqrt {d+e x}+\sqrt {c} (d+e x)\right )}{64 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}}+\frac {3 e \left (2 c^2 d^4+5 a c d^2 e^2+7 a^2 e^4-\sqrt {c} d \sqrt {c d^2+a e^2} \left (2 c d^2+4 a e^2\right )\right ) \log \left (\sqrt {c d^2+a e^2}+\sqrt {2} \sqrt [4]{c} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}} \sqrt {d+e x}+\sqrt {c} (d+e x)\right )}{64 \sqrt {2} a^2 \sqrt [4]{c} \left (c d^2+a e^2\right )^{5/2} \sqrt {\sqrt {c} d+\sqrt {c d^2+a e^2}}} \\ \end{align*}

Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 2.06 (sec) , antiderivative size = 382, normalized size of antiderivative = 0.42 \[ \int \frac {1}{\sqrt {d+e x} \left (a+c x^2\right )^3} \, dx=\frac {\frac {2 \sqrt {a} \sqrt {d+e x} \left (11 a^3 e^3+6 c^3 d^3 x^3+a^2 c e \left (5 d^2+16 d e x+7 e^2 x^2\right )+a c^2 d x \left (10 d^2+d e x+12 e^2 x^2\right )\right )}{\left (c d^2+a e^2\right )^2 \left (a+c x^2\right )^2}+\frac {3 i \left (4 c d^2+10 i \sqrt {a} \sqrt {c} d e-7 a e^2\right ) \arctan \left (\frac {\sqrt {-c d-i \sqrt {a} \sqrt {c} e} \sqrt {d+e x}}{\sqrt {c} d+i \sqrt {a} e}\right )}{\left (\sqrt {c} d+i \sqrt {a} e\right )^2 \sqrt {-c d-i \sqrt {a} \sqrt {c} e}}-\frac {3 i \left (4 c d^2-10 i \sqrt {a} \sqrt {c} d e-7 a e^2\right ) \arctan \left (\frac {\sqrt {-c d+i \sqrt {a} \sqrt {c} e} \sqrt {d+e x}}{\sqrt {c} d-i \sqrt {a} e}\right )}{\left (\sqrt {c} d-i \sqrt {a} e\right )^2 \sqrt {-c d+i \sqrt {a} \sqrt {c} e}}}{32 a^{5/2}} \]

[In]

Integrate[1/(Sqrt[d + e*x]*(a + c*x^2)^3),x]

[Out]

((2*Sqrt[a]*Sqrt[d + e*x]*(11*a^3*e^3 + 6*c^3*d^3*x^3 + a^2*c*e*(5*d^2 + 16*d*e*x + 7*e^2*x^2) + a*c^2*d*x*(10
*d^2 + d*e*x + 12*e^2*x^2)))/((c*d^2 + a*e^2)^2*(a + c*x^2)^2) + ((3*I)*(4*c*d^2 + (10*I)*Sqrt[a]*Sqrt[c]*d*e
- 7*a*e^2)*ArcTan[(Sqrt[-(c*d) - I*Sqrt[a]*Sqrt[c]*e]*Sqrt[d + e*x])/(Sqrt[c]*d + I*Sqrt[a]*e)])/((Sqrt[c]*d +
I*Sqrt[a]*e)^2*Sqrt[-(c*d) - I*Sqrt[a]*Sqrt[c]*e]) - ((3*I)*(4*c*d^2 - (10*I)*Sqrt[a]*Sqrt[c]*d*e - 7*a*e^2)*
ArcTan[(Sqrt[-(c*d) + I*Sqrt[a]*Sqrt[c]*e]*Sqrt[d + e*x])/(Sqrt[c]*d - I*Sqrt[a]*e)])/((Sqrt[c]*d - I*Sqrt[a]*
e)^2*Sqrt[-(c*d) + I*Sqrt[a]*Sqrt[c]*e]))/(32*a^(5/2))

Maple [F(-1)]

Timed out.

\[\int \frac {1}{\left (c \,x^{2}+a \right )^{3} \sqrt {e x +d}}d x\]

[In]

int(1/(c*x^2+a)^3/(e*x+d)^(1/2),x)

[Out]

int(1/(c*x^2+a)^3/(e*x+d)^(1/2),x)

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 5779 vs. \(2 (774) = 1548\).

Time = 3.93 (sec) , antiderivative size = 5779, normalized size of antiderivative = 6.28 \[ \int \frac {1}{\sqrt {d+e x} \left (a+c x^2\right )^3} \, dx=\text {Too large to display} \]

[In]

integrate(1/(c*x^2+a)^3/(e*x+d)^(1/2),x, algorithm="fricas")

[Out]

Too large to include

Sympy [F(-1)]

Timed out. \[ \int \frac {1}{\sqrt {d+e x} \left (a+c x^2\right )^3} \, dx=\text {Timed out} \]

[In]

integrate(1/(c*x**2+a)**3/(e*x+d)**(1/2),x)

[Out]

Timed out

Maxima [F]

\[ \int \frac {1}{\sqrt {d+e x} \left (a+c x^2\right )^3} \, dx=\int { \frac {1}{{\left (c x^{2} + a\right )}^{3} \sqrt {e x + d}} \,d x } \]

[In]

integrate(1/(c*x^2+a)^3/(e*x+d)^(1/2),x, algorithm="maxima")

[Out]

integrate(1/((c*x^2 + a)^3*sqrt(e*x + d)), x)

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1664 vs. \(2 (774) = 1548\).

Time = 0.42 (sec) , antiderivative size = 1664, normalized size of antiderivative = 1.81 \[ \int \frac {1}{\sqrt {d+e x} \left (a+c x^2\right )^3} \, dx=\text {Too large to display} \]

[In]

integrate(1/(c*x^2+a)^3/(e*x+d)^(1/2),x, algorithm="giac")

[Out]

-3/32*(2*(a^2*c^2*d^4*e + 2*a^3*c*d^2*e^3 + a^4*e^5)^2*(sqrt(-a*c)*c*d^3*e + 2*sqrt(-a*c)*a*d*e^3)*abs(c) + (2
*a^2*c^4*d^8*e + 9*a^3*c^3*d^6*e^3 + 19*a^4*c^2*d^4*e^5 + 19*a^5*c*d^2*e^7 + 7*a^6*e^9)*abs(a^2*c^2*d^4*e + 2*
a^3*c*d^2*e^3 + a^4*e^5)*abs(c) + (4*sqrt(-a*c)*a^3*c^6*d^13*e + 25*sqrt(-a*c)*a^4*c^5*d^11*e^3 + 67*sqrt(-a*c
)*a^5*c^4*d^9*e^5 + 98*sqrt(-a*c)*a^6*c^3*d^7*e^7 + 82*sqrt(-a*c)*a^7*c^2*d^5*e^9 + 37*sqrt(-a*c)*a^8*c*d^3*e^
11 + 7*sqrt(-a*c)*a^9*d*e^13)*abs(c))*arctan(sqrt(e*x + d)/sqrt(-(a^2*c^3*d^5 + 2*a^3*c^2*d^3*e^2 + a^4*c*d*e^
4 + sqrt((a^2*c^3*d^5 + 2*a^3*c^2*d^3*e^2 + a^4*c*d*e^4)^2 - (a^2*c^3*d^6 + 3*a^3*c^2*d^4*e^2 + 3*a^4*c*d^2*e^
4 + a^5*e^6)*(a^2*c^3*d^4 + 2*a^3*c^2*d^2*e^2 + a^4*c*e^4)))/(a^2*c^3*d^4 + 2*a^3*c^2*d^2*e^2 + a^4*c*e^4)))/(
(a^4*c^5*d^9 + 4*a^5*c^4*d^7*e^2 + 6*a^6*c^3*d^5*e^4 + 4*a^7*c^2*d^3*e^6 + a^8*c*d*e^8 - sqrt(-a*c)*a^4*c^4*d^
8*e - 4*sqrt(-a*c)*a^5*c^3*d^6*e^3 - 6*sqrt(-a*c)*a^6*c^2*d^4*e^5 - 4*sqrt(-a*c)*a^7*c*d^2*e^7 - sqrt(-a*c)*a^
8*e^9)*sqrt(-c^2*d - sqrt(-a*c)*c*e)*abs(a^2*c^2*d^4*e + 2*a^3*c*d^2*e^3 + a^4*e^5)) + 3/32*(2*(a^2*c^2*d^4*e
+ 2*a^3*c*d^2*e^3 + a^4*e^5)^2*(c^2*d^3*e + 2*a*c*d*e^3)*abs(c) + (2*sqrt(-a*c)*a*c^4*d^8*e + 9*sqrt(-a*c)*a^2
*c^3*d^6*e^3 + 19*sqrt(-a*c)*a^3*c^2*d^4*e^5 + 19*sqrt(-a*c)*a^4*c*d^2*e^7 + 7*sqrt(-a*c)*a^5*e^9)*abs(a^2*c^2
*d^4*e + 2*a^3*c*d^2*e^3 + a^4*e^5)*abs(c) + (4*a^3*c^7*d^13*e + 25*a^4*c^6*d^11*e^3 + 67*a^5*c^5*d^9*e^5 + 98
*a^6*c^4*d^7*e^7 + 82*a^7*c^3*d^5*e^9 + 37*a^8*c^2*d^3*e^11 + 7*a^9*c*d*e^13)*abs(c))*arctan(sqrt(e*x + d)/sqr
t(-(a^2*c^3*d^5 + 2*a^3*c^2*d^3*e^2 + a^4*c*d*e^4 - sqrt((a^2*c^3*d^5 + 2*a^3*c^2*d^3*e^2 + a^4*c*d*e^4)^2 - (
a^2*c^3*d^6 + 3*a^3*c^2*d^4*e^2 + 3*a^4*c*d^2*e^4 + a^5*e^6)*(a^2*c^3*d^4 + 2*a^3*c^2*d^2*e^2 + a^4*c*e^4)))/(
a^2*c^3*d^4 + 2*a^3*c^2*d^2*e^2 + a^4*c*e^4)))/((a^4*c^5*d^8*e + 4*a^5*c^4*d^6*e^3 + 6*a^6*c^3*d^4*e^5 + 4*a^7
*c^2*d^2*e^7 + a^8*c*e^9 - sqrt(-a*c)*a^3*c^5*d^9 - 4*sqrt(-a*c)*a^4*c^4*d^7*e^2 - 6*sqrt(-a*c)*a^5*c^3*d^5*e^
4 - 4*sqrt(-a*c)*a^6*c^2*d^3*e^6 - sqrt(-a*c)*a^7*c*d*e^8)*sqrt(-c^2*d + sqrt(-a*c)*c*e)*abs(a^2*c^2*d^4*e + 2
*a^3*c*d^2*e^3 + a^4*e^5)) + 1/16*(6*(e*x + d)^(7/2)*c^3*d^3*e - 18*(e*x + d)^(5/2)*c^3*d^4*e + 18*(e*x + d)^(
3/2)*c^3*d^5*e - 6*sqrt(e*x + d)*c^3*d^6*e + 12*(e*x + d)^(7/2)*a*c^2*d*e^3 - 35*(e*x + d)^(5/2)*a*c^2*d^2*e^3
+ 44*(e*x + d)^(3/2)*a*c^2*d^3*e^3 - 21*sqrt(e*x + d)*a*c^2*d^4*e^3 + 7*(e*x + d)^(5/2)*a^2*c*e^5 + 2*(e*x +
d)^(3/2)*a^2*c*d*e^5 - 4*sqrt(e*x + d)*a^2*c*d^2*e^5 + 11*sqrt(e*x + d)*a^3*e^7)/((a^2*c^2*d^4 + 2*a^3*c*d^2*e
^2 + a^4*e^4)*((e*x + d)^2*c - 2*(e*x + d)*c*d + c*d^2 + a*e^2)^2)

Mupad [B] (verification not implemented)

Time = 12.64 (sec) , antiderivative size = 9035, normalized size of antiderivative = 9.82 \[ \int \frac {1}{\sqrt {d+e x} \left (a+c x^2\right )^3} \, dx=\text {Too large to display} \]

[In]

int(1/((a + c*x^2)^3*(d + e*x)^(1/2)),x)

[Out]

atan(((((3*(14336*a^9*c^3*e^11 + 4096*a^5*c^7*d^8*e^3 + 18432*a^6*c^6*d^6*e^5 + 38912*a^7*c^5*d^4*e^7 + 38912*
a^8*c^4*d^2*e^9))/(2048*(a^10*e^8 + a^6*c^4*d^8 + 4*a^9*c*d^2*e^6 + 4*a^7*c^3*d^6*e^2 + 6*a^8*c^2*d^4*e^4)) -
((d + e*x)^(1/2)*(-(9*(16*a^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4
+ 210*a^8*c^2*d^3*e^6 - 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 - 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4
096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2
*d^2*e^8)))^(1/2)*(4096*a^9*c^4*d*e^10 + 4096*a^5*c^8*d^9*e^2 + 16384*a^6*c^7*d^7*e^4 + 24576*a^7*c^6*d^5*e^6
+ 16384*a^8*c^5*d^3*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4
)))*(-(9*(16*a^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2
*d^3*e^6 - 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 - 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e
^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(
1/2) - ((d + e*x)^(1/2)*(441*a^4*c^3*e^10 + 144*c^7*d^8*e^2 + 612*a*c^6*d^6*e^4 + 1089*a^2*c^5*d^4*e^6 + 990*a
^3*c^4*d^2*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*
(16*a^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6
- 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 - 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^1
0*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)*1i -
(((3*(14336*a^9*c^3*e^11 + 4096*a^5*c^7*d^8*e^3 + 18432*a^6*c^6*d^6*e^5 + 38912*a^7*c^5*d^4*e^7 + 38912*a^8*c
^4*d^2*e^9))/(2048*(a^10*e^8 + a^6*c^4*d^8 + 4*a^9*c*d^2*e^6 + 4*a^7*c^3*d^6*e^2 + 6*a^8*c^2*d^4*e^4)) + ((d +
e*x)^(1/2)*(-(9*(16*a^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210
*a^8*c^2*d^3*e^6 - 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 - 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(
a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*
e^8)))^(1/2)*(4096*a^9*c^4*d*e^10 + 4096*a^5*c^8*d^9*e^2 + 16384*a^6*c^7*d^7*e^4 + 24576*a^7*c^6*d^5*e^6 + 163
84*a^8*c^5*d^3*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(
-(9*(16*a^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*
e^6 - 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 - 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 +
a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)
+ ((d + e*x)^(1/2)*(441*a^4*c^3*e^10 + 144*c^7*d^8*e^2 + 612*a*c^6*d^6*e^4 + 1089*a^2*c^5*d^4*e^6 + 990*a^3*c^
4*d^2*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*(16*a
^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 - 21*
c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 - 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6
*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)*1i)/((3*(
144*c^6*d^7*e^3 + 684*a*c^5*d^5*e^5 + 882*a^3*c^3*d*e^9 + 1233*a^2*c^4*d^3*e^7))/(1024*(a^10*e^8 + a^6*c^4*d^8
+ 4*a^9*c*d^2*e^6 + 4*a^7*c^3*d^6*e^2 + 6*a^8*c^2*d^4*e^4)) + (((3*(14336*a^9*c^3*e^11 + 4096*a^5*c^7*d^8*e^3
+ 18432*a^6*c^6*d^6*e^5 + 38912*a^7*c^5*d^4*e^7 + 38912*a^8*c^4*d^2*e^9))/(2048*(a^10*e^8 + a^6*c^4*d^8 + 4*a
^9*c*d^2*e^6 + 4*a^7*c^3*d^6*e^2 + 6*a^8*c^2*d^4*e^4)) - ((d + e*x)^(1/2)*(-(9*(16*a^5*c^5*d^9 - 49*a^2*e^9*(-
a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 - 21*c^2*d^4*e^5*(-a^15*c)^(1/2
) + 105*a^9*c*d*e^8 - 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2
+ 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)*(4096*a^9*c^4*d*e^10 + 4096*a^5*c^8
*d^9*e^2 + 16384*a^6*c^7*d^7*e^4 + 24576*a^7*c^6*d^5*e^6 + 16384*a^8*c^5*d^3*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8
+ 4*a^7*c*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*(16*a^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^(1/2)
+ 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 - 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c
*d*e^8 - 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^
4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2) - ((d + e*x)^(1/2)*(441*a^4*c^3*e^10 + 144*c^7*d
^8*e^2 + 612*a*c^6*d^6*e^4 + 1089*a^2*c^5*d^4*e^6 + 990*a^3*c^4*d^2*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c
*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*(16*a^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6
*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 - 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 -
54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4
+ 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2) + (((3*(14336*a^9*c^3*e^11 + 4096*a^5*c^7*d^8*e^3 + 18432
*a^6*c^6*d^6*e^5 + 38912*a^7*c^5*d^4*e^7 + 38912*a^8*c^4*d^2*e^9))/(2048*(a^10*e^8 + a^6*c^4*d^8 + 4*a^9*c*d^2
*e^6 + 4*a^7*c^3*d^6*e^2 + 6*a^8*c^2*d^4*e^4)) + ((d + e*x)^(1/2)*(-(9*(16*a^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^
(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 - 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*
a^9*c*d*e^8 - 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^
12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)*(4096*a^9*c^4*d*e^10 + 4096*a^5*c^8*d^9*e^2
+ 16384*a^6*c^7*d^7*e^4 + 24576*a^7*c^6*d^5*e^6 + 16384*a^8*c^5*d^3*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*
c*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*(16*a^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^
6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 - 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 -
54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^
4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2) + ((d + e*x)^(1/2)*(441*a^4*c^3*e^10 + 144*c^7*d^8*e^2 +
612*a*c^6*d^6*e^4 + 1089*a^2*c^5*d^4*e^6 + 990*a^3*c^4*d^2*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c*d^2*e^6
+ 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*(16*a^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7
*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 - 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 - 54*a*c*d
^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^
13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)))*(-(9*(16*a^5*c^5*d^9 - 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*
d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 - 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 - 54*a*
c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10
*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)*2i + atan(((((3*(14336*a^9*c^3*e^11 + 4096*a^5*c^7*d^8*e^3 + 1
8432*a^6*c^6*d^6*e^5 + 38912*a^7*c^5*d^4*e^7 + 38912*a^8*c^4*d^2*e^9))/(2048*(a^10*e^8 + a^6*c^4*d^8 + 4*a^9*c
*d^2*e^6 + 4*a^7*c^3*d^6*e^2 + 6*a^8*c^2*d^4*e^4)) - ((d + e*x)^(1/2)*(-(9*(16*a^5*c^5*d^9 + 49*a^2*e^9*(-a^15
*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 + 21*c^2*d^4*e^5*(-a^15*c)^(1/2) +
105*a^9*c*d*e^8 + 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 1
0*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)*(4096*a^9*c^4*d*e^10 + 4096*a^5*c^8*d^9
*e^2 + 16384*a^6*c^7*d^7*e^4 + 24576*a^7*c^6*d^5*e^6 + 16384*a^8*c^5*d^3*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*
a^7*c*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*(16*a^5*c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(1/2) + 8
4*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 + 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e
^8 + 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^
6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2) - ((d + e*x)^(1/2)*(441*a^4*c^3*e^10 + 144*c^7*d^8*e
^2 + 612*a*c^6*d^6*e^4 + 1089*a^2*c^5*d^4*e^6 + 990*a^3*c^4*d^2*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c*d^2
*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*(16*a^5*c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4
*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 + 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 + 54*a
*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 1
0*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)*1i - (((3*(14336*a^9*c^3*e^11 + 4096*a^5*c^7*d^8*e^3 + 18432*
a^6*c^6*d^6*e^5 + 38912*a^7*c^5*d^4*e^7 + 38912*a^8*c^4*d^2*e^9))/(2048*(a^10*e^8 + a^6*c^4*d^8 + 4*a^9*c*d^2*
e^6 + 4*a^7*c^3*d^6*e^2 + 6*a^8*c^2*d^4*e^4)) + ((d + e*x)^(1/2)*(-(9*(16*a^5*c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(
1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 + 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a
^9*c*d*e^8 + 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^1
2*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)*(4096*a^9*c^4*d*e^10 + 4096*a^5*c^8*d^9*e^2
+ 16384*a^6*c^7*d^7*e^4 + 24576*a^7*c^6*d^5*e^6 + 16384*a^8*c^5*d^3*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c
*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*(16*a^5*c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6
*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 + 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 +
54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4
+ 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2) + ((d + e*x)^(1/2)*(441*a^4*c^3*e^10 + 144*c^7*d^8*e^2 +
612*a*c^6*d^6*e^4 + 1089*a^2*c^5*d^4*e^6 + 990*a^3*c^4*d^2*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c*d^2*e^6
+ 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*(16*a^5*c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*
e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 + 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 + 54*a*c*d^
2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^1
3*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)*1i)/((3*(144*c^6*d^7*e^3 + 684*a*c^5*d^5*e^5 + 882*a^3*c^3*d*e^9 +
1233*a^2*c^4*d^3*e^7))/(1024*(a^10*e^8 + a^6*c^4*d^8 + 4*a^9*c*d^2*e^6 + 4*a^7*c^3*d^6*e^2 + 6*a^8*c^2*d^4*e^
4)) + (((3*(14336*a^9*c^3*e^11 + 4096*a^5*c^7*d^8*e^3 + 18432*a^6*c^6*d^6*e^5 + 38912*a^7*c^5*d^4*e^7 + 38912*
a^8*c^4*d^2*e^9))/(2048*(a^10*e^8 + a^6*c^4*d^8 + 4*a^9*c*d^2*e^6 + 4*a^7*c^3*d^6*e^2 + 6*a^8*c^2*d^4*e^4)) -
((d + e*x)^(1/2)*(-(9*(16*a^5*c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4
+ 210*a^8*c^2*d^3*e^6 + 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 + 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4
096*(a^15*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2
*d^2*e^8)))^(1/2)*(4096*a^9*c^4*d*e^10 + 4096*a^5*c^8*d^9*e^2 + 16384*a^6*c^7*d^7*e^4 + 24576*a^7*c^6*d^5*e^6
+ 16384*a^8*c^5*d^3*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4
)))*(-(9*(16*a^5*c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2
*d^3*e^6 + 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 + 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e
^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(
1/2) - ((d + e*x)^(1/2)*(441*a^4*c^3*e^10 + 144*c^7*d^8*e^2 + 612*a*c^6*d^6*e^4 + 1089*a^2*c^5*d^4*e^6 + 990*a
^3*c^4*d^2*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*
(16*a^5*c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6
+ 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 + 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^1
0*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2) + ((
(3*(14336*a^9*c^3*e^11 + 4096*a^5*c^7*d^8*e^3 + 18432*a^6*c^6*d^6*e^5 + 38912*a^7*c^5*d^4*e^7 + 38912*a^8*c^4*
d^2*e^9))/(2048*(a^10*e^8 + a^6*c^4*d^8 + 4*a^9*c*d^2*e^6 + 4*a^7*c^3*d^6*e^2 + 6*a^8*c^2*d^4*e^4)) + ((d + e*
x)^(1/2)*(-(9*(16*a^5*c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^
8*c^2*d^3*e^6 + 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 + 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^1
5*c*e^10 + a^10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8
)))^(1/2)*(4096*a^9*c^4*d*e^10 + 4096*a^5*c^8*d^9*e^2 + 16384*a^6*c^7*d^7*e^4 + 24576*a^7*c^6*d^5*e^6 + 16384*
a^8*c^5*d^3*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9
*(16*a^5*c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6
+ 21*c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 + 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^
10*c^6*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2) + (
(d + e*x)^(1/2)*(441*a^4*c^3*e^10 + 144*c^7*d^8*e^2 + 612*a*c^6*d^6*e^4 + 1089*a^2*c^5*d^4*e^6 + 990*a^3*c^4*d
^2*e^8))/(64*(a^8*e^8 + a^4*c^4*d^8 + 4*a^7*c*d^2*e^6 + 4*a^5*c^3*d^6*e^2 + 6*a^6*c^2*d^4*e^4)))*(-(9*(16*a^5*
c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 + 21*c^2
*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 + 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6*d^
10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)))*(-(9*(16*a
^5*c^5*d^9 + 49*a^2*e^9*(-a^15*c)^(1/2) + 84*a^6*c^4*d^7*e^2 + 189*a^7*c^3*d^5*e^4 + 210*a^8*c^2*d^3*e^6 + 21*
c^2*d^4*e^5*(-a^15*c)^(1/2) + 105*a^9*c*d*e^8 + 54*a*c*d^2*e^7*(-a^15*c)^(1/2)))/(4096*(a^15*c*e^10 + a^10*c^6
*d^10 + 5*a^11*c^5*d^8*e^2 + 10*a^12*c^4*d^6*e^4 + 10*a^13*c^3*d^4*e^6 + 5*a^14*c^2*d^2*e^8)))^(1/2)*2i - (((d
+ e*x)^(1/2)*(6*c^2*d^4*e - 11*a^2*e^5 + 15*a*c*d^2*e^3))/(16*a^2*(a*e^2 + c*d^2)) - (e*(d + e*x)^(3/2)*(9*c^
3*d^5 + 22*a*c^2*d^3*e^2 + a^2*c*d*e^4))/(8*a^2*(a*e^2 + c*d^2)^2) - (3*c*e*(c^2*d^3 + 2*a*c*d*e^2)*(d + e*x)^
(7/2))/(8*a^2*(a*e^2 + c*d^2)^2) + (c*e*(d + e*x)^(5/2)*(18*c^2*d^4 - 7*a^2*e^4 + 35*a*c*d^2*e^2))/(16*a^2*(a*
e^2 + c*d^2)^2))/(c^2*(d + e*x)^4 + a^2*e^4 + c^2*d^4 + (6*c^2*d^2 + 2*a*c*e^2)*(d + e*x)^2 - (4*c^2*d^3 + 4*a
*c*d*e^2)*(d + e*x) - 4*c^2*d*(d + e*x)^3 + 2*a*c*d^2*e^2)